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ABSTRACT: Terminal oxidases catalyze four-electron 
reduction of oxygen to water, and the energy harvested 
is utilized to drive the synthesis of adenosine triphosphate. 
While much effort has been made to design a catalyst 
mimicking the function of terminal oxidases, most 
biomimetic catalysts have much lower activity than native 
oxidases. Herein we report a designed oxidase in 
myoglobin with an O, reduction rate (52 s~') comparable 
to that of a native cytochrome (cyt) cbb, oxidase (50 s~’) 
under identical conditions. We achieved this goal by 
engineering more favorable electrostatic interactions 
between a functional oxidase model designed in sperm 
whale myoglobin and its native redox partner, cyt b., 
resulting in a 400-fold electron transfer (ET) rate 
enhancement. Achieving high activity equivalent to that 
of native enzymes in a designed metalloenzyme offers 
deeper insight into the roles of tunable processes such as 
ET in oxidase activity and enzymatic function and may 
extend into applications such as more efficient oxygen 
reduction reaction catalysts for biofuel cells. 


vo normally exhibit much higher activities 
under milder conditions than small inorganic catalysts, 
even though the metalloenzymes use only a limited number of 
physiologically available metal ions or metal-containing 
cofactors. Rational design of metalloenzymes that mimic native 
enzymes allows the elucidation of structural features and 
mechanisms responsible for the high activity and efficient use of 
earth-abundant metal cofactors and facilitates the design of 
more robust and cost-effective catalysts for many applications.’ 
A key measure of success in this endeavor is the activity of the 
designed enzymes compared with the native enzymes that they 
model. While tremendous progress has been made in the field, 
most reported metalloenzymes display activities that are far 
below those of the native enzymes that they mimic. 

Excellent examples of metalloenzymes of interest are 
terminal oxidases in the respiratory chain, in which O, is 
reduced to H,O in a four-electron (4e7) process and the energy 
harvested from this reduction is converted into a proton 
gradient that drives the synthesis of adenosine triphosphate, the 
universal source of energy for most biological processes.” In 
addition to their important role in bioenergetics, terminal 
oxidases are also of interest to biofuel cell research because the 
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oxidases are among the most efficient oxygen reduction 
reaction (ORR) catalysts, having the lowest overpotentials 
and using earth-abundant metal ions such as iron and copper, 
unlike the best industrial ORR fuel cell catalysts, which require 
precious platinum.° To understand the efficiency of this 
important class of metalloenzymes, many synthetic models 
have been made to mimic their structures and functions, with 
particular focus on the heme-copper center at the active site of 
most oxidases.’ Despite years of effort by many groups, the 
activities of the biomimetic compounds in homogeneous 
solution are still far below those of native oxidases." 

As an alternative approach to biomimetic studies of native 
enzymes using small organic molecules, we use small, stable 
proteins such as myoglobin (Mb) as scaffolds to design 
structural features that mimic those of native enzymes.*’ Mb 
functions as an O, storage protein in biology. To transform Mb 
into an oxidase, we have introduced conserved structural 
elements of heme-copper oxidases, one Tyr and two His 
residues, into the distal pocket of Mb. One such variant, L29H/ 
F43H/G65Y Mb (called G6SY-Cu,Mb), was able to reduce O, 
to water with a catalytic oxygen reduction rate of 0.30 s“! 
(Figure S1) and >1000 turnovers using ascorbate as a reductant 
and N,N,N’,N’-tetramethyl-p-phenylenediamine (TMPD) as a 
redox mediator.“ 

While the above results demonstrate that we can rationally 
design a functional enzyme with a high number of turnovers, 
the catalytic rate (0.30 s~') is still far below that of native 
enzymes (e.g., cytochrome (cyt) cbb, oxidase from Rhodobacter 
sphaeroides has an O, reduction rate of 52 s”').'° Our 
preliminary investigations of G6SY-CugMb and similarly 
designed Mb model enzymes suggest that electron transfer 
(ET) into the catalytic center may be rate-limiting. This finding 
is consistent with studies of native cyt c oxidase’’ and some 
synthetic models showing that electron delivery into the active 
site might be a rate-limiting step. ™”? In searching for new 
redox partners of Mb to increase the ET rate and potentially 
the catalytic activity, we were excited to learn of the work by 
Hoffman and co-workers, who replaced three negatively 
charged amino acids (Asp44, Asp60, and Glu85) in myoglobin 
with three positively charged lysines to improve electrostatic 
interactions between the Mb and its redox partner, cyt bs, 
which is negatively charged.'? Such a D44K/D60K/E8S5K 
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variant, called Mb(+6), enhanced the rate of ET from the 
triplet state of Zn-deuteroporphyrin-substituted Mb to ferric 
cyt b; by >2 orders of magnitude (from 5.5 x 10° to ~1.0 x 10° 
s').'** Inspired by this work, we report herein that the D44K/ 
D60K/E85K mutant of G6SY-Cu;Mb (named G6SY-Cu,Mb- 
(+6) hereafter; see Figure la) exhibited dramatically enhanced 
oxidase activity, from 0.30 to 52 s™!, making it a designed 
metalloenzyme with activity comparable to its target native 


enzymes. 
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Figure 1. (a) Structures of G6SY-Cu,;Mb(+6), showing the 
engineered lysines in blue, and cyt b; (PDB IDs 1CYO for cyt bs? 
and 4FWY for F33Y-CugMb;* rendered through VMD°). (b) Oxidase 
activity of G6SY-Cu,;Mb(+6) in comparison with those of native cyt 
cbb, oxidase and G6SY-CugMb at the same concentration under the 
typical conditions of: NADH, 2 mM; cyt b; reductase, 80 nM; cyt bs, 5 
HM; G65Y-Cu,Mb(+6), 50 nM. The black arrows indicate the 
addition of reductant and the double arrow shows the injection of 
native cyt cbb, oxidase. 


The construction, expression, and purification of G6SY- 
CugMb(+6) were carried out using a protocol reported 
previously.’ The mutations were confirmed by DNA 
sequencing, and the corresponding protein containing the 
mutations was corroborated by mass spectrometry (calculated 
MW 17477; observed MW 17476). The UV-vis 
electronic absorption spectra of G6SY-CugMb(+6) in both 
the ferric and ferrous forms are similar to those of G6SY- 
CugMb (Figure S2a) and the EPR spectrum of ferric G6SY- 
Cu,Mb(+6) is similar to that of G6SY-Cu,Mb (Figure S3), 
suggesting that the mutations cause minimal perturbation of the 
overall protein structure and the structure of the heme active 
Site. 

To determine the oxidase activity as measured by the rate of 
O, reduction, we used an O,-selective electrode commonly 
employed to monitor the catalytic activity in native 
oxidases.'”’'* In the presence of 80 nM NADH-cyt b; reductase, 
the physiological redox partner of cyt b," and 5 uM cyt b; the 
concentration of O, in the presence of 50 nM G6SY- 
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CugMb(+6) decreased rapidly (Figure 1b), and the rate of 
O, reduction was 52 s™". The reaction was carried out in 5 mM 
potassium phosphate buffer (pH 6), which has been shown to 
be the optimal pH for O, reduction in this and other CusMb 
variants. The activity decreased at higher ionic strength 
(Figure S4), suggesting that the ionic strength affects the 
electrostatic interactions between Mb and cyt bs, which is 
consistent with the results reported by Hoffman and co-workers 
in the similar system.'° Therefore, low ionic strength (5 mM) 
was used in this study. Similar to G6SY-CugMb,’ which has 
been shown to mimic copper-independent cyt bd oxidase,” the 
addition of copper to G6SY-Cu,Mb(+6) did not increase the 
oxidase activity.’ As a result, no copper was added to the 
reaction mixture. In contrast, G6SY-CugMb without the D44K/ 
D60K/E85K mutations exhibited a rate of 0.3 s™' under the 
same conditions. Furthermore, the protein without the G6SY 
mutation, Cu;Mb(+6), also exhibited a low O, reduction rate 
of 8.1 s`}, which affirms the importance of the G6SY mutation 
for oxidase activity.* Remarkably, the O, reduction rate of 
G65Y-CugMb(+6) (52 s7') is so high that it is similar to that of 
native cbb,; oxidase (50 s™™) measured under the same 
conditions. 

To confirm that the dramatic increase in the oxidase activity 
of G6SY-Cu,Mb(+6) is due to enhanced ET, we first repeated 
the above experiment by systematically removing each protein 
component involved in the O, reduction. As shown in Figure 
SS, minimal O, reduction occurred in the absence of either 
NADH-cyt b; reductase, cyt b, or G6SY-CugMb(+6), 
suggesting that all three proteins are essential for the enzymatic 
activity. Furthermore, we monitored the oxidation of cyt b; in 
the presence of either G6SY-CugMb or G6SY-CugMb(+6) 
through stopped-flow UV-—vis spectroscopy in order to 
estimate the rates of ET between cyt b; and the two Mb 
variants (Figure 2 and Figure S6). We first found that the 
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Figure 2. Studies of ET between cyt b; and G6SY-CugMb or 
G6SYCuBMb(+6). (a) Representative stopped-flow UV—vis spectra of 
cyt b; oxidation catalyzed by G6SY-CugMb(+6) and (inset) time 
traces at characteristic wavelengths. (b) Time traces of cyt b; oxidation 
represented by the absorption at 556 nm, overlaid with global 
spectroscopic fitting (red lines). 
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pseudo-first-order rate constant for ET between 5 uM cyt b; 
and G6SY-CugMb(+6) increased linearly as the concentration 
of G6S5Y-Cug(+6) increased (Figure S7),\* consistent with a 
bimolecular reaction as observed previously. Then we used the 
same concentration of cyt b; as in the oxidase activity assay (5 
HM) but increased the concentration of the Mb variant from 50 
nM to 0.5 uM, which is still a fraction of the cyt b; 
concentration. The pseudo-first-order rate constant for ET 
between cyt b; and G6SY-CugMb(+6) was determined to be 
(1.2 + 0.2) x 10! s', which is ~400-fold faster than that of 
G6SY-CugMb ((3.0 + 0.5) x 107° s+). This 400-fold increase 
is on the same order of magnitude as that reported by Hoffman 
and co-workers (180-fold increase).'* The slight difference is 
attributable to the differences in the Mb mutations and redox 
cofactors (heme vs Zn porphyrin) used in the two studies. 
More importantly, this ET rate enhancement is similar to the 
rate enhancement of O, reduction, which strongly suggests that 
the enhanced ET rate is responsible for the increased oxidase 
activity. The reduction potentials of G6SY-Cug;Mb(+6) and 
G6S5Y-CugMb were determined to be (115 + 11) and (129 + 
5) mV, respectively, by spectroelectrochemical methods 
(Figure S8). Since the reduction potentials of the two mutants 
are similar, these results suggest that the major contributor to 
the increased ET rate is the enhanced electrostatic interaction 
between the negatively charged cyt b; and positively charged 
G6SY-CugMb(+6) rather than a difference in the ET driving 
force. 

To investigate whether O, reduction in the Mb variant is a 
4e process, as in native oxidases, we measured the amount of 
H,O, produced from the reaction using an amperometric H,O, 
sensor. Under the same conditions as in the O, reduction 
experiment shown in Figure 1b, G6SY-CugMb(+6) generated 
H,O, at a rate of 3.5 s™! (Figure 3a), or 7% of the total O, 
reduced, while G6SY-CugMb produced 22% H,O, (Figure 
S1).* Furthermore, since the NADH used in the O, reduction 
experiment is a two-electron reductant, the ratio of the rate of 
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Figure 3. (a) Measurement of O, reduction and H,O, production in 
the O, reduction reaction. (b) Representative kinetic UV—vis spectra 
of reaction solution with 100 „M NADH, sampled at 0.5 s intervals 
over 100 s. The inset shows NADH oxidation and oxygen reduction 
catalyzed by G6SY-CugMb(+6). 
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NADH consumed in the reaction to that of the O, reduced can 
be used to determine whether the O, reduction is a 4e” 
process, as demonstrated in other heme enzymes.'® The 
NADH oxidation was monitored using the absorption at 352 
nm (see Figure 3b).'* The rate of the reaction was compared 
with that of O, reduction using an O,-selective electrode under 
the same conditions as in Figure 3b. As shown in the Figure 3b 
inset, the ratio of the rate of NADH oxidation to that of O, 
reduction was 1.91 + 0.05 for G6SY-Cu,Mb(+6). This result is 
consistent with the low percentage of H,O, directly probed by 
the above electroanalytical method. Together, these two results 
strongly suggest that O, reduction by G6SY-CugMb(+6) is a 
4e process that produces H,O as in native oxidases. Finally, 50 
nM G6S5Y-Cu;Mb(+6) was able to reduce 220 uM O, in the 
reaction system (Figure 1b), corresponding to >4400 turn- 
overs. 

The rational design of proteins that mimic native enzymes 
both structurally and functionally is both exciting and 
challenging because nature has had millions of years to evolve 
enzymes to carry out complex reactions like O, reduction. 
Therefore, any progress in protein design that moves closer to 
the properties of native enzymes is encouraging, with the 
ultimate goal being to match the activities of native enzymes 
while using the models to understand and explore potential 
applications. While most biochemical studies of native enzymes 
use site-directed mutagenesis to find residues that lower the 
activity of the enzyme, biosynthetic modeling seeks to uncover 
residues that result in a gain of function and increase the 
activity; the two approaches therefore complement each other. 
For example, while ET is known to play a critical role in oxidase 
activity, the exact mechanism is not well understood. A report 
by Brunori and co-workers indicated that the internal ET rate 
accounts for the turnover number and that the ET rate and 
turnover number display the same pH and temperature 
dependence.'” Another report by Fabian and co-workers 
provided evidence that the rate-limiting step is the initial 
electron transfer to the catalytic site.” The present results using 
biosynthetic models strongly suggest that the rate of ET into 
the active site plays a critical role in increasing the enzymatic 
activity in a mechanism like that of a native oxidase. A related 
and exciting area of research is to immobilize native enzymes 
and their variants onto electrodes to understand the protein— 
protein and protein—electrode interactions for efficient ET in 
applications such as biofuel cells.'"*' Since our biosynthetic 
models are much smaller and more robust that native oxidases 
and their variants, the ability to design biosynthetic models that 
match the O, reduction activity of native enzymes may open a 
new avenue to explore applications in biofuel cells. 

In conclusion, we have succeeded in engineering Mb, which 
natively binds O, reversibly, into an oxidase that displays the O, 
reduction activity of a native oxidase with a similar rate and 4e7 
reduction process. We achieved the goal by enhancing the ET 
rate through engineering of more favorable electrostatic 
interactions between Mb and its redox partner cyt b, following 
the introduction of similar structural elements (conserved His 
and Tyr) from the native enzyme into Mb. By demonstrating 
our ability to design a protein whose activity matches that of 
the native enzyme, we have obtained a deeper understanding of 
structural features that are important for oxidase activity, which 
may allow engineering of artificial enzymes for biochemical and 
biotechnology applications such as more efficient ORR catalysts 
for biofuel cells. 
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